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Angle-dependent studies of the gap function provide evidence for the coexistence of two distinct 
gaps in hole doped cuprates, where the gap near the nodal direction scales with the superconducting 
transition temperature T c , while that in the antinodal direction scales with the pseudogap temper- 
ature. We present model calculations which show that most of the characteristic features observed 
in the recent angle-resolved photoemission spectroscopy (ARPES) as well as scanning tunneling 
microscopy (STM) two-gap studies are consistent with a scenario in which the pseudogap has a 
non-superconducting origin in a competing phase. Our analysis indicates that, near optimal doping, 
superconductivity can quench the competing order at low temperatures, and that some of the key 
differences observed between the STM and ARPES results can give insight into the superlattice 
symmetry of the competing order. 
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I. INTRODUCTION 

The curious angle-dependence of the gap in the hole 
doped cuprates has been a subject of intense study for 
some time now. Initially, the gap near optimal doping 
was reported to have the ideal | cos (k x a) — cos (k y a) | form 
expected for d— wave superconductivity 1 , and the devi- 
ations observed in underdoped samples were interpreted 
as evidence for the presence of a third-harmonic compo- 
nent in the gap function™. However, it has been found 
recently that the antinodal gap near (it, 0) and the nodal 
gap near (ir/2, ir/2) possess distinctly different doping de- 
pendencies in that the antinodal gap follows the pseudo- 
gap, while the nodal gap scales with T,£&&. Moreover, in 
near-optimally doped Bi 2 Si'2CaCu208 (Bi2212), the gap 
crosses over from being of a pure d— wave form at low 
temperatures to one displaying a pseudogap character 
above T<&£. This has led to further debate as to whether 
the pseudogap is associated with precursor pairin g 8 ' 9 ' 10 ! 11 
or with a competing ordepM^ ' 12 ' 13 ' 14 . Here we explore 
the latter scenario, and show that a model of competing 
order can naturally explain a number of puzzling fea- 
tures of both angle-resolved photoemission spectroscopy 
(ARPES)^££±i^i±&±£ as well as the more recent scan- 
ning tunneling microscopy (STM) 18 ' 19 ' 20 ' 21 experiments. 

We model the pseudogap as a short-range order (SRO) 
which competes with d— wave superconductivity (dSC) 
and possesses the symmetry of the antiferromagnetic 
(AFM) order of the undoped system. Our analysis is 
based mostly on the use of a mean- field t—t'—t" — t'" — U 
Hubbard model of competing AFM and dSC orders, 
and is similar to the one used previously to successfully 
describe a number of properties of the electron-doped 
cuprate o 22 ' 23 . Such a mean field treatment has been 
shown— to mimic AFM short-range order, where the 
Neel temperature Tn approximates the pseudogap on- 
set temperature T* and the AFM gap approximates the 
pseudogap. Notably, a recent study of the optical prop- 
erties of La2- a; Sr : i : Cu04 (LSCO)^ comes to the conclu- 
sion that the cuprates represent the intermediate cou- 



pling case, which would suggest that the cuprates are 
amenable to an approach such as the present one start- 
ing from the weak coupling limit. 



Concerning technical details, we include the supercon- 
ducting (SC) order empirically through a d— wave pairing 
potential V. The staggered magnetisation S at the nest- 
ing vector Q — (tt,tt), which gives the pseudogap US, 
as well as the SC gap A are computed self-consistently 
at all dopings and temperatures. The bare dispersion is 
modeled within a tight-binding approximation using (in 
meV): t = 250, t' = -25, t" = 12, and t'" = 35. These 
values of the hopping parameters are very similar to those 
adduced earlier for LSCO 26 . Values of Hubbard U and 
the pairing potential V have been adjusted to obtain a 
good fit with the experimentally observed Fermi surface 
(FS) arc length and the overall size of the measured SC 
gap2£. Although we focus in this article on the spin den- 
sity wave (SDW) case, we have also investigated other 
ordered phases, including the charge- and d-density wave 
(C/DDW), and the Pomeranchuk mode. 



The paper is organized as follows. Section II discusses 
SDW-based results and the corresponding ARPES data 
on LSCO, while Section III considers the STM data. Sec- 
tion IV compares predictions of the SDW model with 
other candidates for the pseudogap order such as the 
CDW, the DDW and the Pomeranchuk mode. In Sec- 
tion V we comment on the issue of FS arcs vs FS pockets 
in the light of recent quantum oscillation experiments. 
Section VI points out that even though the magnetic 
properties of the cuprates are well-known to be quite 
asymmetric with respect to electron vs hole doping, the 
electronic properties are substantially more electron-hole 
symmetric. A few concluding remarks are presented in 
Section VII. Some of the relevant technical details of our 
modeling are given in Appendices A-C. 
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FIG. 1: (color online) (a) Theoretical and experimental angle 
dependence of various gaps A in underdoped LSCO (x=0.15 ) 
along the FS, where <f> = denotes the antinodal and (f> — 45° 
the nodal direction. Solid lines give our results for the normal 
(T = 50K, red line) and the superconducting state (T = 10K, 
blue line) , and the corresponding experimental data are plot- 
ted with filled symbols of the same coloi 1 ^. Blue dashed line 
shows the computed leading edge gap (LEG) at T = 10K, 
while blue open circles give the corresponding experimental 
LEG as discussed in the text. Blue open triangles denote A 
based on peak positions from the data of Ref. l2a . (b) Com- 
puted energy distribution curves (EDCs) at different momen- 
tum points along the FS (see inset) for the normal (T = 50K, 
blue lines) and the SC state (T = 10K, red lines). Spec- 
tra at the bottom of the figure refer to the antinodal direc- 
tion (<j> = 0), while those at the top to the nodal direction 
(<j> = 45°). Blue and red tick marks on the spectra denote 
total gap values while green lines mark the LEG. (c) Same as 
(b), except that this figure refers to the experimental EDCs 
taken from Figs. 2(a) and (b) of Ref. Htl In order to high- 
light spectral changes, normal state spectra (red) are plotted 
on top of those for the SC state (blue) in several cases, even 
though these pairs of spectra are not taken at exactly the 
same angle <f>. 



II. TWO-GAP SCENARIO AND THE ARPES 
DATA 

ARPES experiments in underdoped LSCO report two 
strikingly different angle dependencies of the gap ; 16 i 28 
but find a natural explanation within our approach. Our 



calculations for underdoped LSCO (x — 0.15) are shown 
in Fig. l(a)-(c). In this case, our self-consistent solution 
yields a staggered magnetization of S = 0.21, which in- 
creases very weakly with T up to T c . The SC gap, i.e. 
A(T = 0) is 13meV with T c = 48K, which is higher 
than the experimental value of 37Ki£, reflecting presum- 
ably the neglect of phase fluctuations 29 in our mean field 
treatment. Interestingly, although the computed ratio 
2A/fcsT c = 6.3 is anomalously large, it is in good agree- 
ment with experiments^. 

We consider first the theoretical results in Figs. 1(a) 
and (b). Fig. 1(a) displays calculated and experimental 
gap values along the Fermi surface as a function of the 
angle </> (inset to Fig. 1(b)), where = corresponds to 
the antinodal direction, and <j> = 45° to the nodal direc- 
tion. Focusing on the energy distribution curves (EDCs) 
of the normal state (red spectra) in (b), the pseudogap 
appears as a broad hump feature with a large gap in 
the antinodal direction (<f) = 0) as marked by red tick 
marks, which decreases to zero by <f> — 18°. For cj> val- 
ues between 18° and 45°, the FS contains an ungapped 
nodal pocket or a FS arc (red curve in (a)). Below T c , 
the EDCs (blue spectra) show an additional sharp peak 
and a leading edge superconducting gap (LEG) over the 
whole FS as marked by green lines, while the hump fea- 
ture (blue tick marks) remains in the antinodal region. 
The presence of a peak-dip-hump feature in the SC state 
clearly reflects the two gap behavior, where the peak fol- 
lows a simple d— wave form while the hump traces the 
pseudogap. Note that all theoretical spectra in (b) have 
been broadened by incorporating the effect of small an- 
gle scattering on the quasiparticles^ 2 -; see Appendix B for 
details. This allows the development of a finite spectral 
weight at Ep and the formation of the leading edge gap, 
even though the underlying quasiparticle states lie well 
below Ep at most momenta^. 

The theoretically predicted gaps derived from the nor- 
mal and SC state spectra of Fig. 1(b) are plotted in Fig. 
1(a), and show good agreement with the corresponding 
experimental data. The characteristics of the evolution 
of the theoretical spectra with FS angle in the presence 
of two gaps, and how these spectra differ between the 
normal and the SC state, as discussed above in connec- 
tion with Fig. 1(b), are also seen in the experimental 
spectral of Fig. 1(c). In particular, our theoretical pre- 
diction that the gap is a pure SC gap up to the tip of 
the FS arc around (p — 18° (see (a)), but that it crosses 
over into becoming a total gap composed of SC and pseu- 
dogap thereafter, gives insight into various experimental 
results reported in the literature- 6,28 -. The peak plotted 
by Shi et. al22- corresponds to our calculated SC peak, 
and as shown in Fig. 1(a) (blue triangles), this gap dis- 
plays a simple d— wave form. In contrast, Terashima et. 
alrS- consider the hump feature and the associated data 
show a two- gap behavior (blue dots in (a)). Nevertheless, 
in the latter datai 6 - which is reproduced in (c), the pres- 
ence of the d— wave LEG can be seen. We have obtained 
values of the LEG from the spectra in (c) and plotted 
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FIG. 2: (color online) Self-consistently computed results for 
LSCO at x = 0.20. 27 (a) Temperature dependence of the SC 
gap (blue line) and the pseudogap (red line) showing reentrant 
behavior discussed in the text, (b) Angle dependence of the 
gap at different temperatures. 
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FIG. 3: (color online) (a)-(c) Spectral intensity at different 
energies u below Ef- Yellow dashed lines mark the magnetic 
zone boundary, (d)-(f) Corresponding simulated q— maps ob- 
tained by convoluting the maps (a)-(c), plotted on a log-scale 
to highlight weak features. Scattering vectors qi and q2, 
shown by arrows, are discussed in the text. Red color de- 
notes the highest and blue the lowest intensity. 



these as open circles in (a). A similar two-gap behavior 
can be seen even more clearly in (Bi,Pb)2(Sr,La)2Cu06+5 
(Bi2201) data in Fig. 3 of Ref.[3 

In the underdoped region, the pseudogap is large com- 
pared to the SC gap. The pseudogap increases weakly 
with T while the SC gap decreases with T. This re sults in 
a nearly constant total gap A = y/ A 2 AFM + A| c , which 
is consistent with experimental observations^. But at 
higher dopings, where the size of the SC gap becomes 
comparable to that of the pseudogap, a more interesting 
temperature dependence can emerge as seen in Fig. 2(a). 
Here a pseudogap is present at high temperature, but 
after the SC gap turns on at 130K, 32 the pseudogap 
is suppressed to zero at 50K and thereafter the total 
gap becomes a pure SC gap. This leads to the evolu- 
tion of the angle dependence of A shown in Fig. 2(b). 
At low temperature (T = 20K, black line), the gap is 
pure d— wave, but at high temperature (T = 155K, red 
line), it is a pure pseudogap, and at intermediate temper- 
atures the gap shows a two-gap behavior similar to that 
of Fig. 1(a). Such a transition from a pure d— wave pair- 
ing gap to a pure pseudogap through a region in which 
both gaps coexist has recently been observed in ARPES 
measurements^ on Bi2212 above the optimal doping re- 
gion. The crossover to pure d— wave form at low tem- 
peratures has been taken as evidence that the pseudogap 
is a precursor SC gap. In contrast, our analysis demon- 
strates that the appearance of a pure d— wave form at 
low temperatures can simply be the replacement of one 
kind of order by a more stable competitor. 

Very different temperature and angle-dependencies in 



the underdoped vs optimal/overdoped regions discussed 
above can be readily understood. The pseudogap, which 
originates here from an ordered phase, only partially gaps 
the FS at any finite doping x. The SC gap, on the other 
hand, opens everywhere, except at the nodal points, so 
that if superconductivity is strong enough it can quench 
a preexisting ordered pseudogap^ 3 -. While it is natural to 
take the competing order to be an SDW in electron-doped 
cuprates, where long-range Neel order persists up to op- 
timal doping, the choice is less clear for hole-doping. Ac- 
cordingly, in Section IV below, we explore other choices 
for the pseudogap. 



III. TWO GAP SCENARIO AND THE STM 
DATA 

We turn next to discuss STM results, where also a 
two-gap behavior has been reported recentl y 18 ! 19 ' 20 ! 21 ' 34 . 
In STM one measures the so-called q— map, from which 
the underlying FS and the angle-dependence of the gap 
can be extracted by interpreting the q— map as the 
Fourier transform of a quasiparticle interference (QPI) 
patter n 18 ' 19 . Figure 3 analyzes the relationship between 
q— maps and the one-particle spectral weights within 
our model, where the spectral weights characterize the 
ARPES spectra to the extent that the ARPES matrix 
element 3 ^ can be neglected. The computations are based 
on the LSCO parameters of Fig. 1 for x — 0.15 at 
T = 10K. In Figs. 3(a)-(c), the computed spectral weight 
A(k, ui) is seen to reside mostly in the momentum region 
of 'bananas' or 'arcs' below the AFM zone boundary^. 
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FIG. 4: (color online) Apparent gap A as a function of the FS 
angle (red dots) as obtained from an analysis of the simulated 
q— maps of Figs. 3(d)-(f). Along the red dashed line the 
peaks in the q— maps are quite broad and may be hard to see 
experimentally. Blue solid line, which gives the total gap, as 
well as the dashed blue line giving the SC gap, is reproduced 
from Fig. 1(a) for reference. Inset: White trace depicts the 
normal state FS on which the FS points obtained from the 
analysis of q— maps of Figs. 3(d)-(f) are shown as red dots. 
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FIG. 5: (color online) Computed angle dependence of the 
gap A in the normal (solid lines) and superconducting state 
(dashed lines) for various competing orders: CDW (red lines); 
DDW (green lines) ; and Pomeranchuk instability for one axis, 
which is taken to be the y axis (blue lines). Black dotted 
line gives the pure d— wave gap in the absence of any other 
competing order. 



At low oj the spectral weight is further concentrated in 
two bright red spots at the 'tips' of each banana, but at 
higher energies the weight spreads out more uniformly 
over the whole FS arc as seen for example in Fig. 3(c). 
The corresponding q— maps, modeled as a convolution of 
the spectral intensity, i.e. 7 q = J2k A(k., u)A(k + q, ui), 
are shown in Figs. 3(d)-(f), and display intense peaks at 
the special q— vectors, which connect the bright spots in 
A(k,oj) 20 . Two such vectors qi and q2 are marked in 
several panels in Fig. 3 as examples^!. It is striking that 
at high energy in (c) , when the bright tips of bananas in 
A(k, oj) have essentially disappeared, the q— map in (f) 
more or less loses its pattern of well defined peaks as the 
intensity spreads out over a wide region. 

The simulated q— maps of Figs. 3(d)-(f) can be used 
to reconstruct the pattern of bright spots in Figs. 3(a)- 
(c), and to thus obtain the FS and the angle-dependent 
gap, as is done commonly in analyzing STM data. The 
resulting gap A is seen from Fig. 4 (red dots) to yield the 
SC gap in accord with the ARPES data (blue solid line), 
but only up to approximately the edge of the FS arc. 
Interestingly, for higher w's, ARPES follows the pseu- 
dogap up to the edge of the Brillouin zone boundary at 
= (blue solid line), but the STM-derived gaps remain 
within the AFM boundary up to the end of the FS arc 
at cf> = 18°. As a result the apparent STM gap adduced 
from the q— map shoots up nearly vertically at <j) = 18° 
following the red dashed line. But with increasing energy 
above the maximum of the SC gap of 13meV at <p = 0, 
it becomes difficult to extract A values as the q— map 
gradually loses its well defined peak pattern except at 
q = 0. Note also that the FS points deduced from the 
q— maps stop near the AFM zone boundary (see inset 
to Fig. 4). These key characteristics of the FS and the 



angle-dependence of the gap in Fig. 4 are in remarkable 
accord with the behavior reported in a recent STM study 
of Bi2212 2 i, and reflect the effect of loss of structure in 
q— maps with increasing energy, which was pointed out 
above in connection with Fig. 3. 



IV. CDW, DDW AND OTHER ORDERS 

Although we have focused on properties of the SDW 
state in this article, we have also carried out calcula- 
tions on a number of other competing electron-hole or- 
dered phases, including the CDW, the DDW, and the 
Pomeranchuk mode; see Appendix C for technical details. 
Fig. 5, which summarizes our key results, shows that the 
CDW (red lines) and DDW (green lines) orders with the 
same reduced AFM Brillouin zone as the SDW, yield a 
gap symmetry very similar to that of Fig. 1 in the nor- 
mal as well as the superconducting state. A doped spin 
liquid model for the pseudogap^ gives similar results. 
However, a Pomeranchuk mod o 39 i 40 does not show a true 
gap, but splits the van Hove singularity (VHS) along x 
and y axes. We can, however, obtain two-gap results sim- 
ilar to those of Fig. 1, along one axis, as shown by the 
blue lines in Fig. 5, suggesting that in a multi-domain 
sample it might be hard to distinguish this behavior from 
the experimental two-gap data. Finally, we have studied 
a linear antiferromagnetic (LAFM) phased, which has a 
one-dimensional ordering vector (tt, 0), as might be seen 
in a stripe phase. However, we find that the resulting two 
gap structure displays a very different symmetry pattern, 
which is not consistent with experiments. 
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FIG. 6: (color online) (a)-(b) Computed Fermi surface maps 
for Na-CCOC based on the AFM model at x = 0.10 and x = 
0.25. Tight binding model parameters used in the calculations 
are: t = 126, *' = -36, t" = 10 meV, with U = 4.5* for 
x = 0.1 (gives S = 0.3) and U = 3.2* for x = 0.25 (gives 
5 = 0). (c)-(d) After Ref. |H, where the FS's obtained in 
Na-CCOC via ARPES from Ref. |H are shown. Red line in 
(c) depicts the FS of YBCO adduced from Shubnikov-de Haas 
experiments of Ref. |42j at similar doping levels. 



V. ARCS VS POCKETS 

An important issue is whether experiments see a well- 
defined pocket or merely a Fermi arc. The recent 
Shubnikov-de Haas (SdH) experiments find oscillations 
in underdoped YBa 2 Cu 3 6 . 5 and YBa 2 Cu 4 8 (YBCO) 
and argue for the presence of a closed pocket of approx- 
imately the expected size 4 ^^. Note that the AFM 
model clearly predicts a full pocket, but the calculated 
spectral weight resembles the Fermi arc with little inten- 
sity on the shadow side due to the effect of AFM coher- 
ence factors [see Figs. 3(a)-(c) and the inset to Fig. 4]. 
Fig. 6 shows results based on the AFM model for another 
hole doped cuprate, Cas-^Na^CuOaCla (Na-CCOC) 43 . 
In the underdoped system at x = 0.10, areas of the 
possible FS pocket in Na-CCOC and YBCO would be 
similars 42 While the observed pocket in YBCO (red line 
in (c)) is somewhat too small 4 ^ to satisfy Luttinger's the- 
orem, it should be remembered that YBCO has a bilayer 
splitting which unlike Bi2212 is not small in the nodal di- 
rection. Hence two nodal FS pockets are expected, where 
the smaller pocket is presumably easier to observe in a 
quantum oscillation experiment. 4 ' Finally, we note that 
ARPES data from underdoped LSCO has recently de- 
tected weak spectral weight on the shadow side of the 
nodal pocket— see Fig. 2 of Ref. [H for doping x — 0.03, 
and 0.07. Interestingly, Kaul et a/. 49 point out that this 
shadow is consistent with a conventional AFM metal, but 
it is not expected in the exotic holon metal phase. 




FIG. 7: (color online) Comparison of selected parameters for 
electron-doped NCCO/PCCO^ (blue lines) vs hole-doped 
LSCO (red lines), (a) Effective interaction U/t; (b) Self- 
consistent magnetisation S; and, (c) Total AFM gap 2US. 



VI. ELECTRON VS HOLE DOPING 

The question of electron-hole symmetry has been a 
topic of interest in cuprate physics for some time. The 
magnetic properties display a strong asymmetry 5 ^ in that 
long-range AFM order persists up to optimal doping with 
electron doping, but it disappears at quite low doping lev- 
els in the hole doped case as it is replaced by the myste- 
rious pseudogap phase. Moreover, nanoscale phase sepa- 
ration, which is prominent in hole-doped cuprates, seems 
to be largely absent in electron doped materials 52 . Some 
of these differences could be understood in terms of how 
the magnetic susceptibility evolves with electron vs hole 
dopingi^iiS 3 . 

The electronic properties, on the other hand, appear to 
show greater symmetry in that superconductivity arises 
near a possible quantum critical point (QCP) close to op- 
timal doping, which is associated with a crossover from 
small to large Fermi surface (FS) 54 ' 55 . In this connec- 
tion, Fig. 7 compares the doping dependencies of selected 
parameters within the AFM model in electron doped 
Nd^Ce^CuOa (NCCO) and Pr^Ce^CuOs (PCCO) 
(blue lines) and hole doped LSCO (red lines). Here, pa- 
rameters for electron doping are taken from our earlier 
worl* 5 ^. The effective U in Fig. 7(a) decreases almost 
linearly with doping in the very underdoped region in a 
very similar manner for both electron and hole doping. 
At higher doping, the smaller U found for hole doping 
could be associated with larger screening resulting from 
the proximity of the VHS. The self-consistently calcu- 
lated value of the staggered magnetization S, however, 
remains higher than that on the electron doped side over 
the entire doping range as seen in Fig. 7(b). Remarkably, 
the total AFM gap 2US (t = 0.326 eV for electron dop- 
ing) displays electron-hole symmetry, although the QCP 
is slightly higher for hole doping. 

The effect of the VHS is reflected in the superconduct- 
ing properties as well. At x = 0.15, for example, the 
self-consistently calculated values of the superconduct- 
ing order parameters are as follows: for electron (hole) 
doping, V = -134(-78.4) meV and A = 2.6(13) meV. 
That is, even though the interaction V is much weaker 
for the hole-doped case, the SC gap is found to be larger, 
leading to a larger ratio 2A/fcsT c of 6.3 for hole doping 
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compared to 4.1 for electron doping. 

In the electron-doped cuprates the QCP involves two 
topological transitions as a function of doping. At low 
doping, the FS consists of electron pockets near T. As 
doping increases, a second hole-like pocket emerges along 
the nodal direction around 15% doping (in NCCO) as 
the magnetic gap decreases and the lower magnetic band 
crosses the Fermi level. The appearance of this pocket 
can be detected in ARPES, Hall effect, and penetra- 
tion depth measurements at essentially the same dop- 
ing level^. At higher doping, the gap collapses and 
the FS evolves into a single large hole- like sheet. A re- 
verse scenario seems to be followed for hole doping: A 
nodal hole-like FS pocket is present at the lowest dop- 
ings (corresponding to the Fermi arc), while at higher 
dopings Hall effect evidence has been found^ for the 
appearance of electron- like pockets in Bi2201, again at 
approximately 15% doping. In our model, these would 
be the T-centered electron pockets associated with the 
upper magnetic band. 



VII. CONCLUSION 

In conclusion, our model calculations show that most 
of the characteristic features observed in ARPES as well 
as STM two-gap studies are consistent with a scenario 
in which the pseudogap has a non-superconducting ori- 
gin in a competing phase. In contrast, a precursor su- 
perconductor model of the pseudogap will have difficul- 
ties explaining why the higher order harmonic content 
of the gap grows with increasing temperature as seen in 
Fig. 2. Our computed spectra not only show the pres- 
ence of a feature which scales with the pseudogap, but 
also display a superconducting low-energy gap leading 
to a peak-dip-hump structure. As doping increases to- 
ward the quantum critical point of the pseudogap, we 
find that a region of reentrant superlattice order could 
reappear in the system. Our analysis highlights electron- 
hole symmetry of electronic properties of the cuprates, 
even though the magnetic properties are well-known to 
be quite asymmetric. 
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APPENDIX A: MODEL FOR THE 
COEXISTENCE OF SDW AND d-SC ORDERS 

Models of competing SDW and SC order have been 
studied for many year o 33 i 56 ' 57 . In our calculations, we 



use a one band tight binding model Hamiltonian where 
antifcrromagnetism is included via a Hubbard U and su- 
perconductivity as in the BCS theory^ 



H 



^ fret cr + U N ct - Cr+c-. - cr, , 



k,k' 



(Al) 



+ Vv r (fc,fc')ct c f - c f - ct , 
■L^, v ' ' k,] -k,l -k',l fc'.T' 

k,k' 

where ct (cg CT ) is the electronic creation (destruction) 



operator with momentum k and spin a (a is the opposite 
spin). The bare particle dispersion with respect to the 
chemical potential Ep is given by 

i% = -2t[c x {a) + c y (a)\ - Wc x (a)c y (a) 

- 2t"[c x (2a) + c y (2a)} 

- At"'[c x {2a)c y {a) + c x (a)c y {2a)} - E F ,(A2) 



with Ci{aa) = cos{akia) and a being the lattice constant. 
t, t' , t" and t'" are TB hopping parameters. Defining the 
Nambu operator in the magnetic Brillouin zone (MBZ) 
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we can write the above Hamiltonian in the MBZ as 
-US A t \ 
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where the order parameters S and represent the stag- 
gered magnetization at nesting vector Q = (jr, it) and the 
superconducting gap, respectively. In the mean field ap- 
proximation, these are defined by 
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E 

k' 



-k',l 



A 9k 



(A6) 



where the d x i_ y i — orbital phase factor is gj: — [c x (a) — 
c y (a)]/2. We diagonalize the Hamiltonian of Eq. I A4I by 
the Bogolyubov method and the corresponding unitary 
matrix can be easily constructed^, 



Ut 



t 


a k u t 




~ a k v t 






a k uZ k 


°Pk4 




a k v t 










~ a k vZ k 





a k v k 

a k u ~k ) 



(A7) 
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The Bogolyubov coefficients are chosen to be 




\ 



1 ± 



E v - 
k 



(A8) 
(A9) 



(q) 2 + (usr. 



where £f = ± £ £+<3 )/2 and £ Qj? = 
The four resulting quasiparticle bands have energies 
±Et, ±EZ, where 



(AfO) 



and 2/ = ± refers to the upper (+) and lower (-) magnetic 
bands (U/LMB). The 4x4 Matsubara Green's function 

can be defined from the Nambu operator as G(k, r— r') = 



T r 1 S r g(r)'W't(T / ) ), whose Fourier transformation gives, 



G(fc, (j, iuj n ) = U£ a {iuj n - iJdiag) 1 U\ g 



(AH) 



where -ffdiag is the diagonalized Hamiltonian containing 
the eigenvalues in order \Et. EZ —Ei, — EZ]. The cor- 

k k k k 

responding 4x4 spectral function is defined in the stan- 
dard form A(k, a, iw n ) = — Im[G(fc, a, ioj n )]/n. We calcu- 
late the self-consistent values of various order parameters 
for each hole doping x, by simultaneously solving the fol- 
lowing set of equations 

/DO J 
—A n (k,a,u + iS)f(«;UA12) 
-nsn ■"TT 



S = 



/CO 7 
— A 12 (fc >ff , W + t<5)/H, (A13) 



and 
A 



<t,uj +i5)f{w). (A14) 



f(u>) = l/(cxp (uj/kBT) + 1) is the Fermi function at 
temperature T where is the Boltzmann constant. 



APPENDIX B: BROADENING DUE TO SMALL 
ANGLE SCATTERING 

In both ARPES and STM calculations we model the 
large spectral broadening to be due to elastic small angle 
scattering of the Cooper pairs (neglecting pair breaking 
effects), in which case the Green's function remains of the 



same form but with renormalized parameters uj — > Co = 
luZ^(lu) and — » Ag = Aj*Zr(to). The renormalization 
factor is taken to be^, 



'■-Al 

k 



(Bl) 



where is the self energy correction due to impurity 

scattering, which is related to the normal state scattering 
rate as22 



Et(u)=sgn(u)(C +C 1 cj p ). 



(B2) 



The power p = 3/2 is assumed to apply for holes, as found 
for electron doped cuprates^. The other parameters Co 
and Ci are found by fitting to the experimental broad- 
ening. For ARPES calculations Co = C\ — 50meV. For 
STM, the uj— dependence is neglected and Co = 30meV. 



APPENDIX C: OTHER MODELS OF 
COMPETING ORDER 

We have also studied various other ordered phases such 
as the CDW, the DDW, the Pomeranchuk instability, and 
the linear antiferromagnetism (LAFM). All these orders 
possess the same super lattice behavior and/or the pseu- 
dospin character near the Van Hove singularity and thus 
are possible candidates for the origin of the pseudogap. 
When the Hamiltonian is solved in the mean field approx- 
imation, the eigenvalues are similar to Eq. lATOl for SDW 
except that E^ is different for each phase. 39,41 Defining 
the interaction as V t with i = CDW, DDW, Pom, LAFM, 
Eq can be written as 



2£ik, (CI) 



where Ri represent the gap in various phase o 39 ' 41 ; in 
Pomeranchuk mode R x — Vp onl g^, in DDW phase R y = 
VddwS£- I n CDW and LAFM phase the gap is similar 
to the SDW phase except that in CDW the gap does 
not depend on the spin orientation of the system and 
thus R z ,a = Vcdw- The LAFM leads to a similar result, 
with R z . a — cVlafM'S', except here the nesting of the 
FS occurs along Q = (ir,0). In each case the interac- 
tion was adjusted to match the gap and the FS arc for 
LSCO at x = 0.15 [see Fig. 5] The resulting values are 
Vddw = 200, Vcdw = 160 meV, whereas in the Pomer- 
anchuk mode, R x simply renormalizes t to have unequal 
values of x and y, t x / y = t ± 1.2 meV. We find that 
LAFM has the wrong pseudogap symmetry to explain 
the experiments. 
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